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ABSTRACT 

We present the energy-dependent power spectral density (PSD) and cross-spectral properties of Mkn 
766 obtained from a six-revolution XMM-Newton observation in 2005. The resulting PSDs, which have 
highest temporal frequency resolution for an AGN PSD to date, show breaks which increase in temporal 
frequency as photon energy increases; break frequencies differ by an average of ~0.4 in the log between 
the softest and hardest bands. The consistency of the 2001 and 2005 observations 5 variability properties, 
namely PSD shapes and the linear rms-flux relation, suggests the 2005 observation is simply a low-flux 
extension of the 2001 observation. The coherence function is measured to be ~0.6-0.9 at temporal 
frequencies below the PSD break, and is lower for relatively larger energy band separation; coherence 
also drops significantly towards zero above the PSD break frequency. Temporal frequency-dependent 
soft-to-hard time lags are detected in this object for the first time: lags increase towards longer time 
scales and as energy separation increases. Cross-spectral properties are the thus consistent with previous 
measurements for Mkn 766 (Vaughan & Fabian 2003) and other accreting black hole systems. The results 
are discussed in the context of several variability models, including those based on inwardly-propagating 
viscosity variations in the accretion disk. 


1. INTRODUCTION 

Seyfert active galactic nuclei (AGNs) and stellar-mass 
black hole X-ray binary systems (XRBs) both exhibit 
rapid, aperiodic X-ray variability that likely originates 
in the innermost regions of these compact accreting ob- 
jects. The dominant X-ray radiation is generally thought 
to be inverse Comptonization of soft seed photons by a 
hot corona (e.g., Shapiro et al. 1976, Sunyaev & Titarchuk 
1980), though the exact geometry is uncertain and numer- 
ous configurations have been invoked (e.g., Zdziarski et 
al. 2003). The X-ray variability can be characterized by 
fluctuation power spectral density functions (PSDs) which 
show the ” red-noise” nature of the variability at relatively 
high temporal frequencies in both AGNs and XRBs. Mod- 
eling of the broadband PSDs of XRBs usually utilizes 
some combination of one, two, or more Lorentzians, which 
tend to dominate in the so-called low/hard energy spectral 
state, plus a broadband noise component characterized as 
(/) oc / _1 , which tends to dominates in the high/soft state 
(e.g., in Cyg X-l; see Axelsson et al. 2005). 

Broadband PSDs have also been constructed for AGNs, 
allowing modeling of the overall PSD shape. Specifically, 
’breaks’ in the power-law PSD models on time scales of 
a few days or less have been identified for roughly a 
dozen AGN (e.g., Uttley, M c Hardy & Papadakis 2002; 
Markowitz et al. 2003; Vaughan, Fabian & Nandra 2003b, 
hereafter VFN03; M c Hardy et al. 2004, 20.05), with power- 


law slopes <-2 above the break and ~-l below it. The 
similarity in Seyferts’ and XRBs’ broadband PSD shapes 
and scaling of PSD break time scale with black hole mass 
(e.g., Markowitz et al. 2003, M c Hardy et al. 2004), as well 
as the presence of a linear relation between flux and abso- 
lute rms variability (e.g., Uttley & M c Hardy 2001, Edelson 
et al. 2002), support the notion of similar X-ray variabil- 
ity mechanisms being present in both classes of accreting 
compact objects. That is, Seyferts and XRBs are simply 
scaled versions of each other in black hole mass and X-ray 
variability time scale 

Further support .for this picture comes from the ob- 
served energy-dependence of the PSD power-law slope 
above the PSD break. Previous studies of both XRBs (e.g., 
Nowak et al. 1999a, Lin et al. 2000) and Seyferts (Nan- 
dra & Papadakis 2001, VFN03, M c Hardy et al. 2004, and 
Markowitz 2005) have suggested that the slope flattens as 
photon energy increases. In addition, Leighly (2004) found 
the structure function of the narrow line Seyfert 1 galaxy 
1H 0707-495 to flatten towards higher photon energies on 
time scales < 3 ksec. No obvious trends in break frequency 
with energy were reported. 

The observed cross-spectral properties, namely the co- 
herence function (Vaughan & Nowak 1997) and time lags, 
are also similar between Seyferts and XRBs. The coher- 
ence is generally seen to be close to unity over a wide range 
of temporal frequencies in both AGNs and XRBs, though 


1 



2 


Mkn 766 X-ray Variability 


usually dropping significantly below unity above the PSD 
break (e.g., Nowak et ah 1999a, 1999b in XRBs; VFN03, 
M c Hardy et al. 2004 and Markowitz 2005 in AGN). Rel- 
atively harder X-rays are seen to generally lag those at 
softer X-ray energies, with the lag between bands increas- 
ing with increasing energy separation. However, lags are 
are also observed to increase with decreasing temporal fre- 
quency (e.g., Miyamoto & Kitamoto 1989, Miyamoto et al. 
1991, Nowak et al 1999a in XRBs), with time lags usually 
<1-10% of the Fourier time scale. In the low-hard state 
of Cyg X-l, the time lags increase iii a step-like fashion 
(e.g., Nowak 2000), which can be explained if the lag at 
each step is associated with individual Lorentzian com- 
ponents in the broadband PSD. Obtaining time lags via 
the cross-spectrum requires extremely high quality data 
with sufficient sampling and number of points. Results 
for AGN so far have yielded time lag spectra with many 
fewer pointer than in XRBs, but also showing the general 
increase in time lags as temporal frequency decreases (Pa- 
padakis, Nandra & Kazanas 2001, VFN03, M c Hardy et al. 
2004, Markowitz 2005, and Arevalo et al. 2006a). 

Mkn 766 is a well-studied Narrow-Line Seyfert 1, objects 
in which large- amplitude rapid X-ray variability is com- 
monly observed. Timing properties of Mkn 766 have been 
studied previously by Vaughan & Fabian (2003; hereafter 
VF03) , who measured the broadband PSD, finding a break 
from steep PSD power-law slope of <-2.5 to ^-1 near 5e-4 
Hz. VF03 also found the slope above the break to flatten 
from ~-2.7 at soft X-rays to ~-2.2 at hard X-rays; how- 
ever, they saw no obvious trends in break frequency with 
energy. VF03 were the first to measure the coherence in 
Mkn 766. In addition, Vaughan et al. (2003a) observed 
the linear rms-flux relation in Mkn 766. 

In this paper we present the timing properties of Mkn 
766 from a 500 ksec long-look with XMM-Newton in 2005. 
The average 0.2-12 keV flux was ~2.5 x 1.0“ 1:1 erg cm" 2 
s'" 1 (e.g., Turner et al. 2006), Jane/ James/Lance: Is 
this value correct? Which paper should be refer- 
enced here, Turner et al. 06 or Miller et al. 06 ? 
just under half that of the 2001 observation, 5.9 x 10" 11 
erg cm" 2 s"" 1 (e.g., Pounds et al. 2003). We use the enor- 
mous wealth of data available from the 2005 long-look to 
construct PSDs with the highest frequency resolution see 
an in AGN PSD so far, yielding tight constraints on PSD 
break frequency. We confirm the rms-flux relation in the 
2005 observation and compare the 2001 and 2005 PSD 
shapes and rms-flux relations. In one of two companion 
timing papers, Papadakis et al. (in prep.; hereafter Pa- 
per II) revisit the broadband PSD shape, using both data 
from the 2005 longJook and from RXTE monitoring data. 
They confirm the break near 5xl0 -4 Hz seen by VF03, 
but use Lorentzians to directly probe similarities with the 
broadband PSD model shapes of XRBs. In the other pa- 
per, Arevalo et al. (in prep.; Paper III) continue the study 
of the energy-dependent variability properties via Fourier- 
resolved spectroscopy. 

This paper is organized as follows: §2 describes the data 
reduction. §3 describes the analysis, including measure- 
ment of the energy-dependent PSD, variability amplitudes 
and autocorrelation functions, as well as stationarity tests. 
§4 describes the cross-spectral properties. The results are 
discussed in §5 in the context of phenomenological vari- 


ability models. Finally, a brief summary is given in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

Mkn 766 was observed by XMM-Newton on 2005 May 23 
- 2005 June 3, over revolutions 999-1004. This paper uses 
data taken with the European Photon Imaging Camera 
(EPIC), which consists of one pn CCD back-illuminated 
array sensitive to 0.15-15 keV photons (Striider et al. 
2001), and two MOS CCD front-illuminated arrays sen- 
sitive to 0.15-12 keV photons (MO SI and MOS2, Turner 
et al. 2001). Data from the pn were taken in Small Win- 
dow mode; data from both MOSes were taken in Large 
Window mode. The medium filter was used for all de- 
tectors. Light curves' were extracted using SAS version 
7.0 using standard extraction procedures. Data were se- 
lected using event patterns 0-4 for the pn and 0-12 for 
the MOSes. Hot, flickering or bad pixels were excluded. 
Source light curves were extracted from circular regions of 
radius 40"; background light curves were extracted from 
circles of identical size, centered ~3 7 away from the source. 
We checked for pile-up; there was no significant pile-up in 
the pn and at most a few percent pile-up in the MOSes. 
The 2-12 keV pn light curve was analyzed for background 
flares; any data where the pn background rate exceeded 
0.2 ct s" 1 were excluded. Most data excluded occurred at 
the beginnings or ends of revolutions (most notably, the 
final 18 ksec of revolution 999 data, the final 17 ksec of 
revolution 1003, and the final 13 ksec of revolution 1004). 
This screening also introduced a 10 ksec gap in the revo- 
lution 1003 light curve. Because such large gaps tend to 
reduce the statistical significance of variability parameters 
derived. over the full duration, the light curve was split 
into two smaller light curves, hereafter designated 10G3A 
and 1003B, for a total of seven virtually uninterrupted 
light curves. All other gaps introduced by background 
screening were small (< a few ksec) and in the analyses 
below, fluxes during gaps were linearly interpolated from 
adjacent points. Data were extracted over the 0.2-12 keV 
(total), 0. 2-0.7 keV (soft), 0.7-2 keV (medium), and 2-12 
keV (hard) bandpasses (also referred to as T, S, M and H 
bands, respectively, below). Light curves from the three 
EPIC cameras were summed, using only data taken when 
all three cameras were in operation simultaneously. Given 
the wide bands used here, the differences between pn and 
MOS response shapes have a negligible effect: the average 
pn and MOS photon energies for each band differed only 
by small amount. In addition, Mkn 766 is a soft-spectrum 
source, and the inclusion of the MOS data compared to 
the pn data alone significantly improved signal/noise in 
the M and H band. Light curves were binned to 60 s; they 
are shown in Figure 1. Table 1 lists mean count rates and 
net exposure times. Fractional variability amplitudes F var 
were calculated for the S, M and H bands for each light 
curve,* using the formulation of Vaughan et al. (2003a); the 
results, listed in Table 1, confirm relatively high levels of 
variability, as expected for a NLSyl. 

3. THE ENERGY DEPENDENCE OF THE HIGH-FREQUENCY 

PSD 

The high-quantity (uninterrupted and highly variable) 
light curves allows us to create a very high quality PSD 
for Mkn 766; in particular, the temporal frequency resolu- 
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tion allows very tight constraints in break frequency, §3.1 
describes how the PSDs were constructed; §3.2 describes 
model fits to the PSD shapes. 

3.1. PSD Construction 

The construction of the PSDs is summarized here; for 
further details on PSD construction, see e.g., Uttley et 
al. 2002, Markowitz et al. 2003, VFN03, and references 
therein. Paper II explores the broadband fit using 2-11 
keV XMM-Newton data and RXTE Proportional Counter 
Array (PCA) monitoring, but since this paper’s aim is 
to explore the PSD behavior across as wide an energy 
range as possible, we do not consider the RXTE data as 
the PCA bandpass does not go below 2 keV. First, for 
each of the seven light curves, the mean was subtracted. 
Periodograms were then calculated for each of the seven 
light curves separately using a Discrete Fourier Transform 
(DFT; e.g., Oppenheim & Shafer 1975). All periodogram 
points from the seven light curves were combined and 
sorted in temporal frequency. Following Papadakis and 
Lawrence (1993), they were binned in temporal frequency 
in groups of 25. The PSD normalization of van der Klis 
(1997) was used. This process yielded a PSD for each 
of the T, S, M and H bands, each with a usable temporal 
frequency range: 3.2 x 10 -5 — 8.3 x 10~ 3 Hz. Power at tem- 
poral frequencies above ~ 3 x 10” 3 Hz were dominated by 
the white noise power due to Poisson noise; this observed 
level of power agreed well with predicted values calculated 
as Pp sn a= 2 (fi -f B)p“ 2 , where fi and B are the net source 
rate and background rate, respectively. The PSDs P(f) 
are plotted in Figures 2a and 3 (/ - P(f ) space and / - 
/xP(/), respectively). Above the break, and below fre- 
quencies where the power due to Poisson Noise dominates, 
the hard band PSD has the highest amplitude, consistent 
with observations for other AGN and XRBs. 

3.2. PSD Model Fits 

We now discuss model fits to the broadband PSD shape. 
Narrow features in the PSD (e.g., features which, in 
the PSDs of XRBs, are usually modeled with narrow 
Lorentzians and may be identified as quasi-periodic oscil- 
lations) are discussed in Paper II. Monte Carlo simulations 
(e.g., Uttley et al. 2002) are frequently employed in PSD 
measurement when it is necessary to account for PSD mea- 
surement distortion effects and to assign errors to poorly- 
sampled PSD points. However, for the present PSDs, we 
expect no aliasing since each of the seven light curves were 
continuously sampled and large gaps were excluded. The 
effects of red noise leak from temporal frequencies lower 
than those sampled here are likely small given the mea- 
sured PSD shape, and at any rate are assumed not to 
vary significantly with photon energy. Finally, because 
each PSD bin contains a sufficient number of unbinned pe- 
riodogram points (25), PSD errors are normal, removing 
the need for Monte Carlo simulations to determine PSD 
errors. We employ least-squares fitting (e.g., Bevington 
1969) to determine the PSD amplitude in all fits below. 
Errors below correspond to Ax 2 == 2.71. 

We first fit an unbroken power-law of the form 
P(f) = A(f/f oTP + Ppsn 

where the normalization A is the PSD amplitude at the 
frequency / 0 and fi is the power law slope; e.g., positive 


values of fi denote red-noise PSD slopes. Pp sn is the con- 
stant level of power due to Poisson noise, calculated to be 
0.18, 0.34, 0.44 and 0.98 Hz'" 1 for the T, S, M and H bands, 
respectively. Pp sn was added to the model as opposed to 
being subtracted from the data, to avoid the possibility 
of obtaining unphysical negative powers. Results for the 
best-fit model are shown in Table 2; residuals to the fits 
are shown in Figure 2b. The fits are quite poor; large 
residuals are apparent and signal the need for a break in 
the PSD model, e.g., as per VF03. 

We next tried fitting broken power-law models. In AGN 
PSDs, it has not always been clear if the breaks are sharp 
or instead follow a slow bend. 

We fit a sharply-broken power-law model of the form 



A(///b)- 7 + Ppsn, 
A(f/f h )-0 + Pp S n, 


f<fb 
/ > fb 


Here, 7 and fi are the power-law slopes below and above 
the break frequency /t>; A is the PSD amplitude at fb. We 
also tested a slowly-broken power-law model, 


P(f) = Am 1 + U/f b)^” 7 ) + Ppsn 

Here, 7 and fi are the power-law slopes far below or above, 
respectively, the rollover frequency fb. A is the PSD nor- 
malization. 

Results for the best-fit models with /3, 7, /b and A free 
are listed in Table 3; residuals are shown in Figures 2c and 
2d. The somewhat high values of yfi for the T, S and H 
bands are due to a few outlier points above ICC 2 / 5 Hz; 
larger temporal frequency bins would alleviate this prob- 
lem but decrease the temporal frequency resolution ofthe 
PSD. The broken power-law fits are significantly better 
(at >99.8% confidence in an F-test) compared to the un- 
broken power-law model fits. Both sharply- and slowly- 
broken model shapes yield similar fits. There is no obvi- 
ous trend in fi with photon energy in either case. In the 
singly-broken fits, 7 tends to increase as photon energy de- 
creases. fb is consistent with increasing as photon energy 
increases. 

We explored if it was significant to keep fi and 7 thawed 
in the fits. We repeated the fits, with fi and 7 frozen at 
the average value of the S, M and H fits. Results are listed 
in Table 3. The sum of S, M and H x 2 /do/ (degrees of 
freedom) increased from 518.6/432 to 526.8/434 for the 
sharply-broken PSDs, and from 521.5/432 to 529.7/438. 
It is this not significant at greater than the ~1 a level to 
thaw fi and 7. With fi and 7 frozen, /b, which has typical 
errors of 0.07 in the log, is seen to increase with increasing 
photon energy (errors do not overlap at the 90% confidence 
level). Refitting the PSDs with /&, beta and 7 frozen yields 
a total (S+M+H) xfj&of of 550.3/441 (sharply-broken) 
and 561.2/441 (slowly-broken). Compared to the model 
fits with just fi and 7 frozen, we find that it is significant 
at 99.97% (sharply-broken) and >99.99% (slowly-broken) 
confidence in an F-test to keep fb thawed in the fits. 

The data are not able to significantly distinguish be- 
tween a sharp break or a slow rollover, and the limited 
dynamic range of the data precludes testing more compli- 
cated PSD model shapes. However, we conclude that for 
all model fits, regardless of whether the power-law slopes 
were fixed or frozen in the fits, the data are always consis- 
tent with fb increasing with photon energy, fb increases 
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by an average of ^0.3 in the log between and S and M 
bands, and ~0.1 in the log between and M and H bands. 
Additionally, considering the lack of obvious trends of (3 
(in both the sharply- and singly-broken fits) or 7 (in the 
slowly-broken fits), the results seem consistent with a sce- 
nario wherein the PSD shape for all bands is identical, but 
the break frequency increases with photon energy. 

Finally, fractional variability amplitudes F var for the S, 
M and H-band light curves, listed in Table 1, show that the 
medium band is almost always the most variable, a result 
identical to VF03 and similar to VFN03 for MCG-6-30-15. 
For light curves associated with red-noise variability pro- 
cesses, F var will be dominated by the longest-term trends; 
the results here are consistent with Figures 2a and 3 where 
one can see that at temporal frequencies below ~ 2 x 10~ 4 
Hz, the medium band PSD has the highest amplitude. 
An jF V ar spectrum for the 2005 observation is presented in 
Miller et al. (2006) who also demonstrate that the variabil- 
ity in the soft and hard bands is diluted by the presence of 
less- variable components, namely the Compton reflection 
spectrum in the hard band and by the soft excess? and 
narrow emission lines in the soft band. 

4. STATION ARITY TESTS 

Compact accreting black holes display “weakly station- 
ary” behavior in that the mean and variance both show 
scatter over time (e.g., the linear rms-flux relation seen in 
Seyferts and XRBs), although both the expectation value 
of the fractional variability amplitude and the underlying 
PSD shape is expected to remain constant in time, e.g, 
over time scales shorter than decades to centuries in AGN. 

4.1. PSD consistency 

We tested for strong non-stationarity in the form of a 
change in the shape of the underlying PSD, between the 
2001 and 2005 PSDs and between revolutions in the 2005 
long-look. For two PSDs observing the same stationary 
process at two different times, the only difference in PSD 
values at a given temporal frequency should be consis- 
tent with the expected scatter. The scatter in the mea- 
sured PSD about the intrinsic, underlying PSD follows a 
X 2 distribution with two degrees of freedom (e.g., Priest- 
ley 1981); Papadakis & Lawrence (1995) thus outlined a 
method to test for non-stationarity, defining a statistic S 
to quantify differences between two PSDs. For a station- 
ary process, S will be distributed with a mean of 0 and a 
variance equal to 1 (see also VFN03). 

We calculated PSDs for each revolution from 999 to 
1003 A using the summed pn+MOS light curves as well as 
for the XMM-Newton observation in 2001 May (revolution 
265; see VF03 for details). However, while the pn camera 
was also operated in Small Window mode in the 2001 ob- 
servation, the MOS cameras were not in Large Window 
mode. The 2001 pn data were reduced and analyzed in a 
manner similar to the 2005 data, yielding T, S, M, and H- 
band light curves 128.5 ksec in duration. Light curves for 
the T, S, M and H bandpasses were binned to 60 seconds 
and truncated to a common duration of 52.2 ksec; PSDs 
were generated by binning periodogram points by a factor 
of 25. The S statistic was calculated by summing over tem- 
poral frequencies <10~ 2,6 Hz to exclude PSD points that 
were dominated by power due to Poisson noise. For all 


PSD pairs and bandpasses, |5| was found to 0.54 or less. 
Extending the duration to 75.8 ksec to compare revolu- 
tions 265 and 999-1002 and extending the duration to 92.3 
ksec to compare revolutions 265 and 1000-2 also yielded 
similarly small values of |£| (<0.5) for all bandpasses. To 
summarize, all values of (S') were less than la for all light 
curve pairs and bandpasses; we therefore have found, no 
evidence for strong non-stationarity in the PSDs of Mkn 
766 at any bandpass from one revolution to the next in 
the 2005 long-look, or between 2001 and 2005. 

4.1.1. The linear RMS- Flux Relation 

As another test for the presence of strong non- 
stationarity between the 2005 and 2001 observations, we 
measured variability amplitudes in the T-band light curve 
over short time scales. Specifically, for a stationary pro- 
cess, we expect a linear relation between average rms vari- 
ability amplitude and flux, or, equivalently, we expected 
< jFvar > to be constant in flux. Vaughan et al. (2003a) al- 
ready established the presence of a linear rms-flux relation 
in the 2001 observation. 

To facilitate comparison to the 2001 pn data, only pn 
light curves from the 2005 observation were used in this 
analysis. Light curves from both observations were binned 
to 60 s. We calculated the flux, F v ar , and the rms ampli- 
tude every 2000 seconds (subtracting off the variance due 
to Poisson noise). Light curves with only 10 points or 
less were discarded. Flux, rms amplitude, and F va r points 
for all six revolutions were combined; data were sorted 
according to flux and binned in groups of 10 to yield aver- 
aged fluxes of flux, rms amplitude cr rms , and F var . Results 
are shown in Figure 4 for both the 2005 long-look (black 
points) and the 2001 long-look (gray points); the top and 
bottom panels show the <a rms >-flux and <F var >-flux re- 
lations, respectively. Both observations are consistent with 
a linear relation between flux and <a rms >, as shown in top 
panel, and are also consistent with <F var > being constant 
in flux, as shown in the bottom panel. Solid lines show 
the best-fit relations, calculated assuming equal weighting 
to all points. Despite the average lower flux in the 2005 
observation, both best-fit lines in each panel of Figure 4 
have similar slopes and intercepts, suggesting that char- 
acteristics of the variability process have not significantly 
changed since 2001, and the 2005 observation is merely a 
low-flux extension of the 2001 observation. The positive 
x-intercepts of the best-fit lines for the <o r rms >-flux rela- 
tion suggest the presence of a constant component, quali- 
tatively corroborating the principle-component analysis of 
Miller et al. (2006). As the x-intercepts for the 2001 and 
2005 relations are approximately equal, we infer that the 
the magnitude of the constant component has not changed 
significantly between the two observations. 

Having established that the variability properties in the 
2001 and 2005 observation are identical, we constructed T, 
S, M and H-band PSDs by combining the 2001 and 2005 pn 
data. The resulting PSD had slightly higher frequency res- 
olution (~%20 in linear space with the same periodogram 
binning), but the level of power due to Poisson noise in the 
M and H bands was about about twice as high compared 
to the PSDs in §3. We repeated the PSD model fits; we 
found results generally consistent with those obtained in 
§3, suggesting that the inclusion of the MOS data did not 
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significantly bias the fitting results. 

5. CROSS-SPECTRAL PROPERTIES 

The cross-spectrum is used to quantify interband cor- 
relations as a function of temporal frequency; analogous 
to the fact that the PSD is the Fourier transform of the 
ACF, the cross-spectrum is the Fourier transform of the 
CCF. The cross-spectrum is a complex number; coherence 
and time lags can be derived from the squared magnitude 
and argument, respectively. Because these quantities in 
compact accreting objects tend to be temporal-frequency 
dependent, they cannot be studied using cross-correlation 
functions alone. Further details of the cross-spectrum and 
its relevance to XRB and AGN observations can be found 
in e.g., Vaughan & Nowak (1997), Nowak et al. (1999a), 
Papadakis, Nandra & Kazanas (2001) and VFN03. 

5.1. The Coherence Function 

The coherence function 7 2 (/) between two time series 
measures the fraction of mean-squared variability of one 
series that can be attributed to other, as a function of 
temporal frequency. If the two time series are related by a 
simple, linear transfer function, then they will have unity 
coherence. The coherence is defined as the magnitude 
squared of the cross spectrum normalized by the product 
of each light curve’s PSD. Here, we use the discrete ver- 
sions of the cross-spectrum and coherence functions, given 
in e.g., §5.2.2 of VFN03. 

We correct the coherence function for the influence of 
Poisson noise by following Eq. 8 of Vaughan & Nowak 
(1997). We assume that the condition of high coher- 
ence and high powers required to use this formulation are 
valid. Specifically, the intrinsic variability power must be 
at least a factor of a few times fym x Pp sn , where m is 
the coherence binning factor (quantified below); this con- 
dition is thus valid for temporal frequencies up to roughly 
1 - 2 x 10" 3 Hz. 

Coherence was measured for each of the seven light 
curves for the S, M and H bands separately, then com- 
bined, sorted in temporal frequency, and binned in groups 
of m= 30 to form the final coherence function for the en- 
tire 2005 data set. As shown in Figure 5, the coherence 
is near unity (>0.8-0. 9) over temporal frequencies lower 
than 'v 6 x 10“ 4 Hz for the S-M and M-H pairs. For S-H, 
the coherence is much lower, ~0.6-0.8, at these temporal 
frequencies. At temporal frequencies above ~ 6 x 10~ 4 
Hz, the coherence significantly drops sharply, reaching 
7 2 (/) < 0.5 by ~ 8 x 10~ 4 Hz. 

This cutoff at high temporal frequencies can be intrin- 
sic, but could also be due to Poisson noise. Monte Carlo 
simulations were used as a check on the accuracy of the 
coherence estimates and to determine if significant devia- 
tions from unity could be attributed to the effect of Pois- 
son noise. For each of the seven observed light curves, 
two light curves were simulated using the same random 
number seed and assuming identical PSD shapes (the av- 
erage of the best-fitting singly-broken PSD model for each 
bandpass pair) to produce two light curves with intrin- 
sic unity coherence. The simulated light curves were re- 
binned and resampled to match the observed light curve, 
and rescaled to match the mean observed count rates in 
the softer and harder bands. Poisson noise was then added 


to the data according to the Poisson distribution, i.e., data 
were randomly deviated using a Gaussian convolved with 
the square root of the observed mean error squared divided 
by the sampling time. One hundred sets of simulated data 
were produced. For each set, the coherence and its uncer- 
tainty was calculated as for the real data. At most tempo- 
ral frequencies below ~ 10“ 3 Hz, the simulated coherence 
was reasonably close to unity, and the uncertainty in the 
coherence was in reasonable agreement with the scatter, 
indicating that the coherence estimation was reasonably 
accurate. The 90% and 95% lower limits to the coherence 
were determined. Any measured coherence point lying be- 
low these lines thus indicates a drop in coherence that is 
intrinsic to the source, and not an artifact of the Pois- 
son noise, at 90% or 95% confidence. For most temporal 
frequencies below ~ 10~ 3 * 5 Hz, the artificial drop in co- 
herence is small. 

The energy-dependence of the high-frequency PSD slope 
could also contribute to a reduction in measured coher- 
ence. The Monte Carlo simulations were repeated, assum- 
ing the best-fit sharply-broken power-law model shapes. 
The effects of Poisson noise were again added. The simu- 
lated coherence was typically only 1-2%, 2-5% and 1-3% 
lower for the S-M, S-H and M-H bands, respectively, com- 
pared to simulations assuming identical PSD shapes. The 
effect of having differing PSD shapes is thus small. The 
90% and 95% lower limits to the coherence were deter- 
mined and are plotted in Figure 5. For most temporal 
frequencies below ~ 10“ 3,5 Hz, the artificial drop in co- 
herence is again small. However, the observed deviations 
from unity are generally larger than the artificial drop due 
to Poisson noise. The decrease in coherence as the en- 
ergy separation of the bands increases and as temporal 
frequency increases past ~ 6 x 10“ 4 Hz is thus likely in- 
trinsic. 

Having quantified the decrease in coherence e.g., due to 
Poisson Noise at each temporal frequency, we added that 
amount (e.g., as determined by the Monte Carlo simu- 
lations) to obtain an estimate of the ’intrinsic’ coherence. 
Fitting the ’intrinsic’ coherence functions with exponential 
decay functions of the form 7 2 (/) = Cexp(— f/f c ) (/ c is 
the temporal frequency cutoff and C is a constant) ignor- 
ing points above 1.2 xl0“ 4 Hz, and using least-squares 
fitting, we find that log(/ c ) for the S-M, S-H and M-H 
coherence functions are —2.90 ± 0.07, —3.07 ± 0.07 and 
—2.90 ± 0.07, respectively (errors correspond to 5x 2 — 
2.71). 

5.2. Temporal Frequency-Dependent Time Lags 

Cross-correlation functions (CCFs) for XRBs and AGN 
generally show a peak near zero lag but are often asym- 
metrically skewed towards hard lags (e.g., McHardy et al. 
2004), suggesting there is a soft-to-hard time lag of some 
sort, but it is not a simple delay. We measured the CCFs 
for each of the seven summed light curves of the 2005 ob- 
servation as well as for the 2001 pn light curve, Using a 300 
s binning. Because large gaps such as those between revo- 
lutions complicate analysis, we calculated CCFs for each of 
the seven light curves separately. CCFs were measured us- 
ing both the Discrete Correlation Function (DCF; Edelson 
& Krolik 1988) and the Interpolated Correlation Function 
(ICF; White & Peterson 1994); the bootstrap method of 
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Peterson et al. (1998) was used to assign errors to lags. All 
peak lags are consistent with zero, with conservative for- 
mal errors of ±300 sec for revolutions 265 and 999-1003 A, 
±1200 sec for revolution 1003B and ±4500 sec for revo- 
lution 1004. The CCFs also exhibit a slight soft-to-hard 
asymmetric skew. The CCFs thus agree qualitatively with 
those measured for other AGN and XRBs, and the logical 
direction is to explore temporal frequency-dependent time 
lags via the cross-spectrum. 

Obtaining sensible measurements of time lags via the 
cross-spectrum is difficult for observations of AGN, how- 
ever; it requires the combination of a large quantity of 
uninterrupted data with a high sampling time resolution 
and a high variability-to-Poisson noise ratio. VF03 did 
not measure lags for the 2001 data; we also attempted to 
derive time lags myself from the 2001 pn data and did not 
get sensible results. 

We used Eqs. 28 and 29 in §5.3.1 of VFN03 to calculate 
the phase lags <j>(f) at discrete temporal frequencies / from 
the summed light curves in the 2005 observation, binned 
to 60 s; time lags were then obtained as r (/) = </>(/) /2rcf . 
Time lags were calculated for each of the seven light curves 
separately, then combined, sorted in frequency, and binned 
in groups of 50. Errors on phase lags were calculated us- 
ing Eq. 16 of Nowak et al. (1999a); in this formulation we 
used the coherence function appropriate for each binned 
point. Sensible lags were obtained for the 4 lowest tem- 
poral frequency bins, spanning 5 x 10" 4 — 4 x 10“ 3 Hz. 
The results are shown in Figure 6. Errors are large, pri- 
marily because coherence in this temporal frequency range 
is not optimally close to unity. Using Eqn. 16 in Nowak 
et al. (1999a) , we estimated the effective sensitivity limit 
on time lag detections due to Poisson noise, using the 
PSD shape and normalization from best-fit M-band singly- 
broken PSD model shape in §3. Limits .are shown as dotted 
lines in Figure 6. 

Monte Carlo simulations were performed to verify that 
the time lags were not an artifact of Poisson noise nor of 
the differing PSD slopes between two bands. The same 
Monte Carlo simulations run for the coherence function in 
the previous subsection were used; e.g., light curves were 
simulated with zero intrinsic time lag between them. Time 
lags were calculated and binned in a manner identical to 
the real data. For all simulations, the simulated time lags 
were close to zero (deviations of the average time lag from 
zero were <32 s for the lowest temporal frequency bin and 
<13 s for the next three bins). The average error was also 
close to the standard deviation of the estimates. We can 
conclude that the observed time lags are not an artifact of 
Poisson noise or PSD shape differences, and are henceforth 
assumed intrinsic to the source. 

The time lag spectrum displays the general qualities 
seen in many previous time lag spectra: a general increase 
towards lower temporal frequencies and as the energy sep- 
aration of the bands increases. Fitting a relation of the 
form r(f) == Cf a yields values for the power-law slope a 
of -0.46 ±0.18, -0.93 ± 0.17, and -1.26 ±0.62 for S-M, 
S-H and M-H, respectively. Alternatively, assuming that 
the time lags are proportional to r(/) = Cf~ l for sim- 
plicity allows us to calculate fractional time lags, i.e,, the 
lags 5 fraction of the temporal frequency. We find that C == 
0.013 ± 0.003, 0.020 ±0.004, and 0.010 ±0.006 for S-M, S- 


H and M-H, respectively, consistent with closure between 
the three bands. The physical implications of these time 
lags and comparisons to other objects will be discussed 
further in section 5. 


6. DISCUSSION 

An XMM-Newton observation of Mkn 766 spanning six 
revolutions in 2005 caught the source in a relatively low 
flux state; with the average 0.2-12 keV flux just under 
half that for a 128 ksec XMM-Newton observation in 2001. 
The data reduction yielded seven long, virtually uninter- 
rupted light curves of this highly- variable source totalling 
455 ksec of exposure. The resulting PSD has the high- 
est temporal frequency resolution for an AGN PSD to 
date, allowing a breakthrough in constraining the break 
frequency with photon energy, as discussed in §6.1. In ad- 
dition, time lags have detected in Mkn 766 for the first 
time; §6.2 discusses the cross-spectral results and com- 
pares them to those of other compact accreting objects. 
Finally, the results are discussed in the context of phe- 
nomological variability models in §6.3, 

6.1. Summary of PSD Results 

The PSD analysis showed that while breaks are signifi- 
cantly detected at all bands, model fits cannot distinguish 
between a sharp break or a slow rollover despite the high 
temporal frequency resolution. In general, there are no 
obvious trends in the power-law slope below or above the 
break, suggesting that the data could be consistent with a 
PSD shape which is the same at all bands. Importantly, 
however, the high temporal frequency resolution allows us 
to place tight constraints on the break frequency - typi- 
cally ~0.07 in the log, and much tighter constraints than 
VF03 or even VFN03 for the three XMM-Newton revo- 
lutions of data for MCG-6-30-15. The break frequency 
tends to increases by only ~0.4 in the log between the S 
and H bands (a factor of ~10 in average photon energy), 
suggesting that PSD measurement spanning a wider range 
of photon energies would be needed to more easily sepa- 
rate energy-dependent break frequencies. It is also possible 
that this phenomenon may not have been widely reported 
in XRBs because all XRB PSDs rely on RXTE s PC A, 
which does not detect <2 keV photons. 

Most previous Seyfert PSDs simply did not have the 
lower temporal frequency resolution to detect significant 
changes in fb with photon energy such as those seen here. 
It is plausible that for a limited dynamic range PSD, the 
effect of having an underlying PSD whose break frequency 
increasing with photon energy could lead to apparent mea- 
surements of /?, the power-law slope above the break, flat- 
tening with energy. Alternatively, it could be the case 
that this phenomenon might exist only in Mkn 766; high- 
quality, high-resolution PSDs are needed for many other 
Seyferts to determine whether this is the case. 

Finally, we find no evidence for the underlying variabil- 
ity process to have changed since the 2001 observation, 
despite the lower average flux. The PSD shapes between 
2001 and 2005 are consistent and both observations appear 
to have consistent linear rms-flux relations and constant 
<F var > values. 
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6.2, Summary of cross- spectral results and comparison to 

other objects 

The coherence is generally flat at temporal frequencies 
below ~ 6 x 10" 4 Hz, with values typically 0.8 — 0.9 
for S-M and M-H and - 0.6 - 0.8 for S-H. The coher- 
ence drops towards zero at higher temporal frequencies; 
Monte Carlo simulations show these deviations from unity 
to he intrinsic to the source and not the result of Poisson 
noise or energy-dependent PSD shapes. The decrease in 
coherence as energy separation increases has been reported 
previously for Mkn 766 by VF03; we in fact see values of 
the coherence similar to VF03. This similarity again sup- 
ports the conclusion of §3.1, that there is no evidence for 
the underlying variability process to have changed signif- 
icantly between 2001 and 2005 despite the overall lower 
average flux in 2005. 

In other black hole systems, the coherence at temporal 
frequencies below the break is generally between 0.8 and 
1.0, though coherence between bands similar to the S and 
H bands used here can be as low as ~0.6-0.7 (e.g., VFN03 
& McHardy et al. 2004); Mkn 766’s coherence functions 
thus agree with those in other compact objects. 

As explained in Vaughan & Nowak (1997), coherence 
between two energy bands can lost if the transfer function 
relating the light curves is non-linear, or if there are mul- 
tiple, uncorrelated transfer functions. If multiple flaring 
regions exist each contribute to the signal in both bands, 
then it is possible for coherence to be low, even if indi- 
vidual regions yield perfectly coherent variability. This is 
because each region independently requires the same lin- 
ear soft-to-hard transfer function in order to yield unity 
coherence between the final light curves. The fact that co- 
herence in Mkn 766 and several other compact objects is 
not exactly unity may support some degree of non-linearity 
present in the transfer functions relating relatively softer 
and harder bands, and which increases for larger energy 
separations. 

As explained by e.g., Nowak & Vaughan (1996) and 
Nowak et al. (1999a), the sharp drop in coherence at tem- 
poral frequencies above the PSD break could indicate that 
the corona is dynamic on such short time scales. Such 
short time scales could be associated with formation time 
scales or changes in temperature and/or the physical struc- 
ture of the corona. The relatively high levels of coherence 
at low temporal frequencies could indicate that the corona 
is effectively static on long time scales. 

We measure temporal-frequency dependent phase lags 
in Mkn 766 for the first time; lags increase towards lower 
temp frequency and as energy separation of bands in- 
creases, as seen in other black hole systems. Fractional 
time lags are about 0.01-0.02 the Fourier time scale. These 
values are consistent with most values measured so far for 
other AGN (see e.g., Figure 8 of Arevalo et al. 2006a), with 
fractional time lags usually measured to be < 0.2, though 
Ark 564’s time lag spectrum (Arevalo et al. 2006a) seems 
to be an exception. 

Assuming a PSD break frequency of 10” 3 ' 4 Hz, the Mkn 
766 lag spectrum is probing temporal frequencies 0.13-1.00 
times the break frequency. We can compare Mkn 766’s 
fractional lags to those measured for Cyg X-l within ± a 
decade of temporal frequency of the PSD breaks. In the 
low /hard and high/soft state of Cyg X-l, fractional time 


lags tend towards 0.03-0.1 (e.g., Nowak et al. 1999a), while 
in the intermediate state, fractional time lags are higher, 
~0. 1-0.2 (e.g., Pottschmidt et al. 2000). On the basis of 
these fractional lag values, Mkn 766 is likely not a su- 
per massive analog of Cyg X-l in the intermediate-state, 
though we cannot determine on the basis of lags alone 
if analogy with the low/hard or high/soft state is appro- 
priate. Furthermore, the time lag spectrum is not of high 
enough quality to distinguish between an intrinsic lag spec- 
trum which is step-like, as in the intermediate state of Cyg 
X-l, or adheres to a form more resembling a power-law. 

6.3. Phenomenological Variability Models 

Many models have been invoked to explain the red 
noise variability properties of Seyferts and XRBs, includ- 
ing shot-noise models (e.g., Merloni & Fabian 2001), ro- 
tating hot-spots on the surface of the accretion disk (Bao 
& Abramowicz 1996), self-organized criticality (“pulse 
avalanche” models; Mineshige, Ouchi & Nishimori 1994), 
magnetohydrodynamical instabilities in the disk (Hawley 
& Krolik 2001), and inwardly-propagating fluctuations in 
the local accretion rate (Lyubarskii 1997). As noted by Ut- 
tley, McHardy <fc Vaughan (2005), shot-noise models have 
difficulty reproducing the linear rms-flux relation, and as 
noted by Vaughan & Nowak (1997) and §5.2, shot-noise 
models lead to low observed coherence unless each sep- 
arate emitting region independently has the same linear 
transfer function between energy bands. 

Hording Sz Falcke (2004) note that a pivoting power-law 
model seems qualitatively capable of reproducing many 
of the observed variability characteristics in XRBs (and 
by extension, Seyferts). This model yields high coherence 
over a broad range of time scales, ACF peaks which get 
narrower as energy increases, and phase lags which exhibit 
a dependency on temporal frequency in a manner similar 
to those observed. 

Thermal Comptonization of soft seed photons by a hot 
corona is a likely explanation for the X-ray emission. It is 
supported by X-ray energy spectra, and the soft-to-hard 
time lags could be attributed to the time scale for seed 
photons to diffuse through the corona and undergo mul- 
tiple up-scatterings. However, the simplest Comptoniza- 
tion models do not predict time lags that depend on tem- 
poral frequency, as observed. In addition, these models 
predict that relatively harder energy bands, having under- 
gone more scatterings, will have rapid variability washed 
out, e.g., variability amplitudes at time scales shorter than 
the PSD break time scale should decrease with increas- 
ing energy, contrary to observations. However, the en- 
ergy dependence of the PSD could be obtained with a 
corona whose temperature increased towards smaller radii 
(Kazanas, Hua & Titarchuk 1997). 

A model in which inwardly-propagating variations in 
the local mass accretion rate m are responsible for the 
observed X-ray variability seems to be able to explain 
many of the observational results in Seyferts and XRBs 
(Lyubarskii 1997, Churazov, Gilfanov & Revnivtsev 2001), 
Kotov et al. 2001, Arevalo & Uttley 2006; see also VFN03 
and McHardy et al. 2004). In this model, variations at 
a given radius are associated with the local viscous time 
scale, so relatively smaller radii are associated with rela- 
tively more rapid variations. Perturbations in m propa- 
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gate inward, until they reach, and modify the emission of, 
the central X-ray emitting region. The resulting PSD has 
a 1// form across a broad range of temporal frequencies, 
the net sum of variations from a wide range of radii. How- 
ever, if the X-ray emitting region is extended, it essentially 
acts like a low-pass filter on these 1/f variations: with rel- 
atively higher- frequency variations imprinted on the light 
curve originating only from the smallest radii, this pro- 
duces a cutoff in the high-frequency PSD, as is commonly 
observed. Kotov et al. (2001) suggested that the spec- 
trum can be a function of radius, with relatively harder 
emission emanating from smaller radii, and therefore as- 
sociated with more rapid variability. This class of models 
thus predicts several additional properties which are ob- 
served in Seyferts and XRBs, namely the presence of the 
linear rms-flux relation, the dependence of coherence and 
time lags on temporal frequency and on the energy sepa- 
ration of the bands, CCF lags, and short-term variability 
amplitudes increasing with photon energy. 

Arevalo & Uttley (2006) show that if the emissivity in- 
dices differ for each band (e.g., harder X-ray emission is 
more centrally concentrated, as per Kotov et al. 2001), 
and if the PSD break is at least partly due to filtering pro- 
duced by having an extended X-ray emission region, then 
the PSDs will have break frequencies which increase with 
increasing photon energy. We note that the PSDs plotted 
in Figure 3 for Mkn 766 bear a striking qualitative resem- 
blance to the simulated PSDs in Figure 2 of Arevalo & Utt- 
ley (2006) . The bending in the Mkn 766 PSDs is slow, and 
rollover frequencies are relatively close together, so it is not 
straightforward to relate break frequencies to emissivity in- 
dex. However, we note that the ’soft’ and ’hard 5 PSDs in 
Figure 2 of Arevalo & Uttley (2006), simulated assuming 
emissivity indices which differ by 2, yield break frequen- 
cies which are separated by about a decade. In Mkn 766, 
the S- and H-band break frequencies are not as separated, 
suggesting that, in the context of this model, the emis- 
sivity indices for the S- and H- bands differ by less than 
2. Patricia - you mentioned in your 23 June email 
that Phil has an XSPEC routine to fit PSDs and 
lag spectra by varying emissivity indices — do you 
think it’s needed here? Is this something you can 
do easily? Also, is this next statement accurate? 
Similarly, the relatively small fractional time lags observed 
in Mkn 766 (0.01-0.02 for the S-H bands) could also be 
an indication that the difference in emissivity indices for 
the bandpasses used here is relatively small. Patricia - in 
Ark 564, with fractional time lags being so. high, is 
this an indication of large difference in emissivity 
indices in Ark 564? Are energy-dependent PSD 
break frequencies observed in Akn 564? 

(I got the impression from VF03 that they though co- 
herence, especially S— H, was a bit lower compared to other 
objects. I’m not sure I agree. If it is true, though, does 
that mean that damping could play more of a role in Mkn 
766’s disk compared to other objects?) 

7. CONCLUSIONS 

We have analyzed the energy-dependent variability 
properties of the Narrow Line Seyfert 1 Mkn 766 from 


a 455 ksec XMM-Newton observation. The long, uninter- 
rupted, variable light curves allow us to construct a PSD 
with the highest temporal frequency resolution for an AGN 
PSD to date. We find PSD breaks at nearly the same tem- 
poral frequencies as VF03, but the break frequencies /b are 
seen to increase with photon energy, typically by ~0.3 be- 
tween the S and M bands and by ~0.1 between the M and 
H bands. The data cannot distinguish between a sharply 
broken or slowly-bending PSD model shape. There is no 
obvious trend in PSD power-law slopes, and the PSDs are 
consistent with a universal PSD shape whose break fre- 
quency increases with photon energy. Previous Seyfert 
PSDs lacked the temporal frequency resolution to cleanly 
detect the energy dependene in /t>; previous reports of 
high-frequency power-law slope flattening with increasing 
photon energy can again be understood in the context of a 
universal PSD shape whose break frequency increases with 
photon energy. 

Variability properties between the 2005 and 2001 obser- 
vations are consistent with each other, despite a 0.2-12 
keV flux in 2005 which is less than half that of the 2001 
observation. Comparing the 2001 and 2005 PSDs shows 
the PSDs to be consistent within the expected scatter. A 
linear rms-flux relation, or equivalently, constant <i ? var >, 
was reported for the 2001 observation by Vaughan et al. 
(2003a). These relations are confirmed for the 2005 obser- 
vation, and shown to be consistent with the corresponding 
relations from the 2001 observation. 

Cross-spectral properties are qualitatively consistent 
with those measured for black hole systems previously. 
The coherence is generally flat at temporal frequencies be- 
low 'v 6 x 10" 4 Hz, with values typically ~ 0.8 — 0.9 for 
S-M and M-H and ~ 0.6 — 0.8 for S-H. The coherence 
drops towards zero at higher temporal frequencies; Monte 
Carlo simulations indicate that the drop is intrinsic to the 
source and not an artifact of Poisson noise or differing PSD 
shapes. 

Temporal frequency-dependent time lags are measured 
from the cross-spectrum for the first time in Mkn 766. As 
is the case with time lags in other objects, lags increase 
as temporal frequency increase and as energy separation 
of the light curves increases. Fractional time lags are typ- 
ically <0.2, consistent with many other Seyferts as well as 
with the high/soft and low/hard states of Cyg X-l. 

The results were discussed in the context of several vari- 
ability models, including models incorporating inwardly- 
propagating fluctuations in the local mass accretion rate, 
e.g., Lyubarskii (1997). Notably, the observed increase 
in the PSD break frequency with photon energy is qual- 
itatively consistent with the prediction put forth by e.g., 
Arevalo & Uttley (2006), who model variations in a disk 
with relatively harder X-ray emission being more centrally 
concentrated. 


This work has made use of observations obtained with 
XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA member states 
and the US (NASA). 
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0.2-12 keV EPIC light curves, summed over the pn and MOS cameras, binned to 60 s, for the 2005 observations. All panels have 
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Table 1 


Source and derived variability parameters 


Rev. 

Exposure 

(ksec) 

Band 

Mean 
ct s” 1 

Bkgd 
ct s" 1 

Fv ar 
. . (%) 

999-1004 

455.0 

T 

13.44 

0.06 

40.52 ± 0.04 



5 

6.58 

0.03 

37.56 ± 0.06 



M 

4.89 

0.01 

44.15 ± 0.07 



H 

2.22 

0.01 

32.42 ± 0.10 

999 

75.8 

T 

5.65 

0.04 

24.90 ± 0.15 



S 

3.16 

0.02 

18.91 ± 0.20 



M 

1.76 

0.01 

23.54 ± 0.25 



H 

1.30 

0.01 

22.38 ± 0.29 

1000 

93.7 

T 

11.31 

0.07 

34.36 i 0.09 



S 

5.59 

0.03 

30.01 ± 0.13 



M 

4.12 

0.01 

38.37 ± 0.15 



H 

1.96 

0.02 

27.65 ± 0.22 

1001 

93.0 

T 

15.06 

0.06 

21.19 ± 0.08 



S 

7.35 

0.03 

22.28 ± 0.12 



M 

5.50 

0.01 

21.71 ±0.14 



H 

2.35 

0.01 

17.47 ± 0,21 

1002 

92.3 

T 

18.65 

0.06 

18.13 ± 0.07 



S 

9.00 

0.03 

18.59 ± 0.11 



M 

6.85 

0.01 

19.01 ± 0.12 



H 

2.91 

0.01 

17.66 ± 0.19 

1003A 

52.2 

T 

18.04 

0.06 

13.67 ± 0,10 



S 

8.46 

0.03 

13.82 ± 0,15 



M 

6.86 

0.01 

14.28 ± 0.17 



H 

2.81 

0.01 

13.38 ± 0.27 

1003B 

18.6 

T 

10.77 

0.05 

16.60 ± 0.22 



S 

5.33 

0.03 

14.28 ± 0.32 



M 

3.83 

0.01 

18.45 ± 0.37 



H 

1.84 

0.01 

15.45 ± 0.52 

1004 

19.5 

T 

12.13 

0.06 

17.82 ± 0.21 



S 

5.86 

0,03 

12.87 ± 0.31 



M 

4.45 

0.01 

22.29 ± 0.35 



H 

2.04 

0.02 

20.28 ± 0.51 


Note. — Col. (2) lists the next exposure after screening, as described in §2. T, S, M, and H refer to the 0.2-12, 0.2-0.74, 
0.74-2.0, and 2.0-12 keV bandpasses, respectively. Col. (4) lists the mean count rate after background subtraction, summed 
from all three EPIC cameras. Col. (5) is the background rate summed over all three cameras. 


Table 2 

Unbroken Power-law Model Fits to the PSD 


Band 

X 2 /dof 

P 

A 

T 

234.8/146 

o • qo + 0,.08 
^^_ 0.06 

1.03±0.04 

S 

192.7/146 

cy O9+0-10 

0.89±0.03 

M 

209.9/146 

2.20±0.08 

1.13±0.03 

H 

259.2/146 

i.o»_ 0 06 

1.04±0.03 


Note. — A is the PSD amplitude at 10~ 3 * 50 Hz. 
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Table 3 


Sharply- and Slowly Broken Power-law Model Fits to the PSDs 


Band 

X 2 /dof 

7 

i°g(/b) 

0 

A 

Sharply-Broken Power-Law Model Fits 

T 

170.2/144 

1.56i0.16 

-3.38±0.06 

2.84I-- . 

0.98-L0.04 

S 

M 

H 

171.0/144 

142.6/144 

205.0/144 

1.72db0.18 
1.40i0. 16 
1.28i0. 14 

-3.48i0.ll 

-3.35i0.05 

-3.27i0.06 

2.66+°; 1 if 

2.98±0.20 

2.82d:0.28 

1.02±0.03 

1.06±0.03 

0.87zh0.03 

T 

S 

M 

H 

170.6/146 

174.5/146 

143.9/146 

208.4/146 

1.47 (fixed) 
1.47 (fixed) 
1.47 (fixed) 
1.47 (fixed) 

-3.41±0.07 

-3.53±0.08 

-3.36±0.05 

-3.23±0.06 

2.82 (fixed) 
2.82 (fixed) 
2.82 (fixed) 
2.82 (fixed) 

1.06±0.04 

1.20±0.03 

1.05±0.03 

0.75±0.03 

Slowly-Broken Power-Law Model Fits 

T 

S 

173.5/144 

169.1/144 

A * iO - 0.16 

0.80i0.40 

-3.49±0.07 

-3.92i0.10 

SOO-aio 

2.62±0.10 

-2.60dz0 04 
-0.97T0.03 

M 

147.2/144 

1 99 + 0.14 
A * ■ — 0.16 

-3.29±0.06 

3.46±0.26 

-2.87i0.03 

H 

205.2/144 

1.14±0.12 

-3.17±0.07 

^ 40 + 0-40 
— 0.34 

-2.74i0.03 

T 

S 

M 

H 

175.5/146 

176.1/146 

148.0/146 

205.6/146 

1.12 (fixed) 
1.12 (fixed) 
1.12 (fixed) 
1.12 (fixed) 

-3.43±0.D6 

-3.55±0.G7 

~3.38i0.06 

-3.21dz0.07 

3.19 (fixed) 
3.19 (fixed) 
3.19 (fixed) 
3.19 (fixed) 

-2.49i0.04 

-2.49i0.03 

-2.45i0.03 

-2.65i0.03 


Note. — 



Fig. 2 .— a) Power density spectra for the 2005 observation; black, red, green and blue points denote T, S, M and H bands, respectively. 
See text for details on PSD construction. The solid lines show the best-fit slowly-broken power-law model, though the constant level of power 
due to Poisson noise has been omitted from the model for clarity. The PSD fits are consistent with break frequency increasing with photon 
energy, b) Residuals to the best-fit unbroken power-law model, c) Residuals to the best-fit sharly-broken power-law model, d) Residuals to 
the best-fit slowly-bending power-law model (power due to Poisson noise was included in the models in panels b)-d). 
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Mkn 766 X-ray Variability 



Fig. 3. — PSDs for the 2005 observation, plotted in / x P(F) space. 



FlG. 4. — The 2005 pn data (black circles) and the 2001 pn data (gray open squares) are both consistent with a linear relation between rms 
variability amplitude (top) and with constant average fractional variability amplitude <F VbX >. The black solid and gray dashed lines are the 
best-fit relations for the 2005 and 2001 data, respectively. 
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Temp. Freq. (Hz) Temp. Freq. (Hz) Temp. Freq. (Hz) 


Fig. 5. — Coherence functions for S-M, S— H and M-H light curve pairs. See text for details on calculation of the coherence function. The 
thin solid line denotes unity coherence. The dashed and dotted lines denote, respectively, the 90% and 95% confidence limits for spurious 
lack of coherence due to the combination of Poisson noise and the effect of differing PSD shapes. 



Temporal Frequency (Hz) 

Fig. 6. — Time lag spectra for Mkn 766. Positive lags denote soft band variations leading those in the hard band. The dashed lines denotes 
the sensitivity limit for time lag detection; see text for details. The solid lines denote the best-fit relations to r(f) = Cf~ 1 . 





